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1. Introduction 

The supersymmetry technique is a powerful method in random matrix theory and 
disordered systems. For a long time it was thought to be applicable for Gaussian 
probability densities only [U El [3J H]. Due to universality on the local scale of the 
mean level spacing [HI El [7J [8], this restriction was not a limitation for calculating 
in quantum chaos and disordered systems. Indeed, results of Gaussian ensembles are 
identical for large matrix dimension with other invariant matrix ensembles on this scale. 
In the Wigner-Dyson theory [§] and its corrections for systems with diffusive dynamics 
[TO] . Gaussian ensembles are sufficient. Furthermore, universality was found on large 
scale, too [UJ. This is of paramount importance when investigating matrix models in 
high-energy physics. 

There are, however, situations in which one can not simply resort to Gaussian 
random matrix ensembles. The level densities in high-energy physics [TJ] and finance 
[T3] are needed for non-Gaussian ensembles. But these one-point functions strongly 
depend on the matrix ensemble. Other examples are bound-trace and fixed-trace 
ensembles [H], which are both norm-dependent ensembles [15] . as well as ensembles 
derived from a non-extensive entropy principle [T6l [T71 ITS] . In all these cases one is 
interested in the non-universal behavior on special scales. 

Recently, the supersymmetry method was extended to general rotation invariant 
probability densities [THJ [TJJ1 EH EI] • There are two approaches. The first one is the 
generalized Hubbard-Stratonovich transformation [T5], [2T] . With help of a proper Dirac- 
distribution in superspace an integral over rectangular supermatrices was mapped to a 
supermatrix integral with non-compact domain in the Fermion-Fermion block. The 
second approach is the superbosonization formula [HI [20] mapping the same integral 
over rectangular matrices as before to a supermatrix integral with compact domain in 
the Fermion-Fermion block. 

In this work, we prove the equivalence of the generalized Hubbard-Stratonovich 
transformation with the superbosonization formula. The proof is based on integral 
identities between supersymmetric Wishart-matrices and quadratic supermatrices. The 
orthogonal, unitary and unitary-symplectic classes are dealt with in a unifying way. 

The article is organized as follows. In Sec. El we give a motivation and introduce 
our notation. In Sec. El we define rectangular supermatrices and the supersymmetric 
version of Wishart-matrices built up by supervectors. We also give a helpful corollary for 
the case of arbitrary matrix dimension discussed in Sec. [71 In Sees. [Hand El we present 
and further generalize the superbosonization formula and the generalized Hubbard- 
Stratonovich transformation, respectively. The theorem stating the equivalence of both 
approaches is given in Sec. El including a clarification of their mutual connection. In 
Sec. [7J, we extend both theorems given in Sees. H] and [5] to arbitrary matrix dimension. 
Details of the proofs are given in the appendices. 
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2. Ratios of characteristic polynomials 

We employ the notation defined in Refs. [221 [21] . Herm (/3, N) is either the set of N x N 
real symmetric (j3 = 1), NxN hermitian (f3 = 2) or 2N x 2N self-dual (f3 = 4) matrices, 
according to the Dyson-index (3. We use the complex representation of the quaternionic 
numbers EL Also, we define 

H^!i 2 ' 4} • H^- 1,2} 

and 7 = 7172 . 

The central objects in many applications of supersymmetry are averages over ratios 
of characteristic polynomials 



k 2 

n det(#- {E n2 -ie)t 12N ) 
Z klk2 (E-) = J P(H)^ d[H] 

Herm (/3,7V) H det (H - (E nl - ie)t^ N ) 

n=l 

J P(iJ)Sdet- 1 ^ (H ® t^ {kl+k2) - t 12N ® E-) d[H] (2.2) 

Herm (/3,7V) 

where P is a sufficiently integrable probability density on the matrix set Herm (/3, iV) 
invariant under the group 

f O(iV) , = 1 
U (/3) (iV) = I XJ(N) , (3 = 2 . (2.3) 
[ USp (2JV) , /3 = 4 

Here, we assume that P is analytic in its real independent variables. We use the 
same measure for d[H] as in Ref. [22] which is the product over all real independent 
differentials, see also Eq. (14.111) . Also, we define E = diag (En, . . . , E kl i, E\ 2 , . . . , E k2 2)® 
1^ and E~ = E — lel^^+k-z)- 

The generating function of the fc-point correlation function flUi [27J [T3J, [21] 

p fc (x) = 7 2 _fc / p(^) n tr 6 ( x p - h ^ ri \ ( 2 - 4 ) 

Herm (/3,7V) P=1 

is one application and can be computed starting from the matrix Green function and 
Eq. (I2.2p with k\ = k 2 = k. Another example is the n-th moment of the characteristic 
polynomial [281 ESI EZ] 

Z n (x,(JL)= J P(H)Q(H)det n (H-Et r2k )d[H} } (2.5) 

Herm (f3,N) 

where the Heavyside-function for matrices 0(H) is unity if H is positive definite and 
zero otherwise. [21] 

With help of Gaussian integrals, we get an integral expression for the determinants 
in Eq. (12.21) . Let Aj be the Grassmann space of j-forms. We consider a complex 
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272 Nk 2 

Grassmann algebra [30] A = Aj with 72iV/c 2 pairs {Cjn, C* n }> 1 — n — ^2, 1 < J ' < 

72^, of Grassmann variables and use the conventions of Ref. [22] for integrations over 
Grassmann variables. Due to the Z^-grading, A is a direct sum of the set of commuting 
variables A and of anticommuting variables A 1 . The body of an element in A lies in A 
while the Grassmann generators are elements in Ai. 

Let 1 be the imaginary unit. We take pairs {zj n , z* n }, 1 < n < ki, 1 < j < 

72^, of complex numbers and find for Eq. (12. 2p 



Z klk2 (E-) = (2tt 



n 2 N{k 2 -k{) 72 Nk 2 



TP{K) exp (-zStr BE~) d[(]d[< 



(2.6) 



ki 72iV k 2 72 JV 

where d[z\ = ]J fl dz Jp dz* p , d[Q = [] fl (<Kjpd£j P ) and £ = x A 272 vfc2- The 

p=lj=X p=l j=l 

characteristic function appearing in (12.61) is defined as 



TP{K) 



P(H)exp{itx HK) d[H]. 



(2.7) 



Herm (f3,N) 



The two matrices 
K = 



and B = -VV ] (2.8) 

7 7 

are crucial for the duality between ordinary and superspace. While K is a 7 2 A^ x 72A" 
ordinary matrix whose entries have nilpotent parts, B is a 7(^1 + k 2 ) x 7(^1 + /C2) 
supermatrix. They are composed of the rectangular 72A" x j(ki + /c 2 ) supermatrix 

0|^ 2 = (z u ..., z kl ,Yz* u . . . , F4 i5 Ci, • • • , Ck 2 , YQ, • • • , YCJ, 
V\ m = (zl ...,z* kl ,Y Zl ,..., Yz kl , -C, . . . , -e k2 ,Yb, . . . , Y( k2 f, 

V^\(3 =2 = (*i,...,* fcl ,Cl,...,C*a)> 

^ =2 =(z 1 *,...,z^ 1 ,-Cr,...,-C 2 ) T (2.9) 

The transposition "T" is the ordinary transposition and is not the supersymmetric 
one. However, the adjoint "f" is the complex conjugation with the supersymmetric 
transposition "Ts" 



a 



T S 







T S 








0~12 




T 
°-ll 


T 
°21 


0~21 


0~22 




T 

_ ~°~12 


T 
°~22 _ 



(2.10) 



where a is an arbitrary rectangular supermatrix. We introduce the constant 72^ x 7 2 A" 
matrix 



Y 




1 

4 



Y 



1 
-1 



The crucial duality relation [TBJ [21] 

trK m = StiB m , m G N, 



(2.11) 



(2.12) 
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holds, connecting invariants in ordinary and superspace. As TP inherits the rotation 
invariance of P, the duality relation (I2.12p yields 



Zk lk2 (E-) = (2tt 



I l 2 N(k 2 -k 1 ) l ~ l2 Nk 2 



$(5) exp (-zStr BE~) d[£\d[z 



(2.13) 



Here, $ is a supersymmetric extension of a representation TP® of the characteristic 
function, 

= ^P (StrS m |m G N) = .FP (tr K m \m E N) = FP(K). (2.14) 

The representation TPq is not unique [3]]. However, the integral (12.131) is independent 
of a particular choice [2"Tj . 

The supermatrix B fulfills the symmetry 

YBYl ,/3e{l,4}, 
YB*Y T ,0 = 2 



B* 



(2.15) 



with the supermatrices 








lfc! 







y s <g> 1^ 












Ifc! 








, *V=4 = 










(2.16) 










Y s ® l fe2 







n 2 








1/3=2 



and y 

^1/3=1 = 

and U 1/3=2 



g(t) 

/3,7l fc l:7 2 fc 2 



where 



1 

7! 



lfe 1+ jt 2 and is self-adjoint for every (3. Using the 7r/4-rotations 
t kl t kl 

-ilk! «ljfci , U\[)_ I 

l 2fc2 



1 

7! 
















lfc 2 




(2.17) 







«lfc 2 





lA, 1 +fc 2 ) ^ = UBU^ lies in the well-known symmetric superspaces [32] 



o- G Mat(7/ci/7fc 2 ) 



(7 



(7. 



(7 



yyi ,72 fc 2 ^y^ i72 /c 2 

1 k\k 2 I k k 
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/3=1 



1 
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o y 
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13=2 



p+q 



and y 



/3=4 



y 



o 



i P o 
1„ 



(2.18) 



.(2.19) 



The set Mat(p/g) is the set of (p + q) x (p + g) supermatrices on the complex Grassmann 

2pq 

algebra ® Aj. The entries of the diagonal blocks of an element in Mat(p/g) lie in A 

j=0 

whereas the entries of the off-diagonal block are elements in A 1 . 

The rectangular supermatrix = V'U' is composed of real, complex or 
quaternionic supervectors whose adjoints form the rows. They are given by 



9\ 



({V2Rez jn , VZlmZjn}^^ , {( jn , (* n } 

{Cjn} 

Kn<k 2 j 



l<n<k 2 



Zjn 
Z j+N,n 



— Z 



j+N,n 



Z 



J" 



l<n<ki 



A-) A+) 
A-)* A+)* 



1, 

2. 



(2.20) 



l<n<fc 2 / 
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respectively, where = « ( - 1±1 ^ 2 (Cjn ± Q+N,n)/^- Then, the supermatrix B acquires 
the form 

~ 1 N 

S =~E^4 ( 2 - 21 ) 

' 3=1 

The integrand in Eq. (I2.13P 

F (b \ — $ exp f-zStr EB\ (2.22) 
comprises a symmetry breaking term, 

3 U E U (/3) (7ifci/7 2 fc 2 ) that fffl)/J? {iJBU^j , (2.23) 
according to the supergroup 

f UOSp (+ >(2Jfei/2Jfe 2 ) , /3 = 1 
U (/3) (7i^i/72^)= < U(fci/A; 2 ) , = 2 . (2.24) 

[ UOSp ( - } (2^/2^) , /3 = 4 

We use the notation of Refs. [Ml [22] for the representations UOSp ^ of the supergroup 
UOSp . These representations are related to the classification of Riemannian symmetric 
superspaces by Zirnbauer [32]. The index "+" in Eq. (12.241) refers to real entries in the 
Boson-Boson block and to quaternionic entries in the Fermion-Fermion block and "— " 
indicates the other way around. 



3. Supersymmetric Wishart— matrices and some of their properties 

We generalize the integrand (I2.22p to arbitrary sufficiently integrable superfunctions on 
rectangular (72C+71G?) x (7 2 a+7i&) supermatrices V on the complex Grassmann-algebra 

2(ad+bc) 

A = Aj. Such a supermatrix 
3=0 

v = (*<?, ■ ■ ■ , *£\*S\ • • • = (*g>- • ■ ■ , , *g>*, • • • (3.1) 

is defined by its columns 
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or by its rows 



^(R)t 



{ x nj}x<n<a ' {Cy"> ^ n i}i<„<fe) 
Kn<a ' Kn<6 / 



"nj2 



2>nj'2 %nj 1 
— V* . 



Kn<a 




"J 
Cry 



^nj'2 



Kn<b/ 



Kn<bS 



<n<b. 



{-Xnj,Xnj} l<n<a ,{ynj} 

Kn<6 / 



= 

P 



1, 
2. 



(3.4) 



2. 

4 



(3.5) 



which are real, complex and quaternionic supervectors. We use the complex Grassmann 
variables Xmn and ( mn and the real numbers x mn and y mn . Also, we introduce the 
complex numbers z mn , z mn , z mn i and z mn \. The (7 2 c + 71 d) x (72C + 71 <f) supermatrix 
P = ^~ l VV^ can be written in the columns of as in Eq. (I2.2ip . As this supermatrix 
has a form similar to the ordinary Wishart-matrices, we refer to it as supersymmetric 
Wishart-matrix. The rectangular supermatrix above fulfills the property 



V* = Y rA VY 



ab ' 



(3.6) 

The corresponding generating function (12.21) is an integral over a rotation invariant 
superfunction P on a superspace, which is sufficiently convergent and analytic in its 
real independent variables, 



where 

E~ = diag (E u ® 1 



P( ( j)Sdet- 1 ^ [a® II 



f(c) 

L 2il> 



1 



720+716 



E-)d[a], (3.7) 



?(-vo 



72 1 



T et ? 0{t) 
ijet ^/3,a6 

opt v(-'Z') 



^(t) 



cl <8> I72) ^12 



Si 



■-71 ) 



(12 



1 71 ) — «elL 



72C+7ia!- 



(3. 



be a subset of ^' ab . The entries of elements in lie in A and Ai. The 
fl^E^Iln^ is the Wick-rotated set of by the generalized Wick- 



rotation n 



(R) 



e l ^' 2 lL 71 f,). As in Ref. [21J, we introduce such a rotation for 



f(c) 

L 2ip 



diag (1 72C , e^l 71( f) is also a 



1/) /3,a6 —i)> 

q - diag (1 720 , 
the convergence of the integral (13. 7p . The matrix II 
Wick-rotation. 

In the rest of our work, we restrict the calculations to a class of superfunctions. 
These superfunctions has a Wick-rotation such that the integrals are convergent. We 
have not explicitly analysed the class of such functions. However, this class is very large 
and sufficient for physical interests. We consider the probability distribution 

P(a) = /( ( x)exp(-Stra 2m ), (3.9) 

where m G N and / is a superfunction which does not increase so fast as exp(Strcr 2m ) 
in the infinty, in particular 



limP (ee lc V) = lim exp (-ee"*Str a 2m ) = 



(3.10) 
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for every angle a G [0, 2-71"]. Then, a Wick-rotation exists for P. 

To guarantee the convergence of the integrals below, let % = nf ] V, = V^uf ] 
and = BH\ . Considering a function / on the set of supersymmetric Wishart- 
matrices, we give a lemma and a corollary which are of equal importance for the 
superbosonization formula and the generalized Hubbard-Startonovich transformation. 
For 6 = 0, the lemma presents the duality relation between the ordinary and superspace 
( I2.12j) which is crucial for the calculation of ( 12. 2ft . This lemma was proven in Ref. [20] 
by representation theory. Here, we only state it. 

Lemma 3.1 

Let f be a superfunction on rectangular supermatrices of the form Ii3. 1\) and invariant 
under 

f(%,Vlj=f(v i ,U\UV^) , (3.11) 



for all V and U G \J {f3 \a/b). Then, there is a superfunction F on the \J^\c/d)- 
symmetric supermatrices with 

F{%) = f{%,V^) . (3.12) 

The U ^(c/cQ-symmetric supermatrices are elements of S^ afc . The invariance 
condition (13.111) implies that / only depends on the rows of by tyffi^ms for arbitrary 
n, m, r and s. These scalar products are the entries of the supermatrix V^V}_^ which 
leads to the statement. 

The corollary below is an application of integral theorems by Wegner [31] worked 
out in Refs. [351 136] and of the Theorems III. 1, III. 2 and III. 3 in Ref. [22]. It states that 
an integration over supersymmetric Wishart-matrices can be reduced to integrations 
over supersymmetric Wishart-matrices comprising a lower dimensional rectangular 
supermatrix. In particular for the generating function, it reflects the equivalence of 
the integral ( 13. 7\} with an integration over smaller supermatrices [22] . We assume that 
a = a-2(b-b)/[3>0 with 

P = 4 and b G 2N + 1 
else 

Corollary 3.2 

Let F be the superfunction of Lemma V3.1\ real analytic in its real independent entries 
and a Schwartz- function. Then, we find 

J F{%)d[V\ = C J F{B^)d[V] (3.14) 

where B = j^VV . The sets are VI = R^+^d/p x A 2(ad+bc) a nd *R = mP mU ^ x 
^-2(ad+bc)> the constant is 



c 



" 7i" 


(b-b)c 


"72" 


(a—a)d 


. 2 . 




. 2 . 





(3.15) 
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and the measure 
d[V\ 



1 [ dXmnl j | dy mni J [ dCrandCmn J [ ^mn^X™ ■ ( 3 - 16 ) 



l<m<a 
Kn<c 



l<m<b 
Kn<d 



l<m<fe 
Kn<c 



l<m<a 
Kn<<2 



l<Z</3 l<K4//3 

The (72C+71C?) x (72S+716) supermatrix V and its measure d[V] is defined analogous to V 
and d[V], respectively. Here, x mna and y mna are the independent real components of the 
real, complex and quaternionic numbers of the supervectors and \IQ 2 , respectively. 

Proof: 



We integrate F over all supervectors and \E^ J except . Then 



,(R) 



,(R) 



(3.17) 



only depends on tfJJ^JJ*. The integration set is 9t' = ^a(c-i)+4bd/ /3 x A 2{ad+b{c 



and the measure d[V^n] is d[V] without the measure for the supervector ty^'. With 
help of the Theorems in Ref. [Ml 1351 1361 122] . the integration over is up to a 
constant equivalent to an integration over a supervector . This supervector is equal 
to * n in the first a-th entries and else zero. We repeat this procedure for all other 
supervectors reminding that we only need the invariance under the supergroup action 
U (/3) (b-b/b-Pj on / as in Eq. f[3~TO embedded in U (/3) (a/6). This invariance is 
preserved in each step due to the zero entries in the new supervectors. □ 
This corollary allows us to restrict our calculation on supermatrices with 6=1 only 
to /3 — 4 and 6 = for all f3. Only the latter case is of physical interest. Thus, we give 
the computation for b = in the following sections and consider the case b = 1 in Sec. 
[3 For b = we omit the Wick-rotation for B as it is done in Refs. [151 ETJ due to the 
convergence of the integral (13.71) . 



-1)) 



4. The superbosonization formula 

We need for the following theorem the definition of the sets 
<7i rj rf 



Si 



pq 



-2,pq 



J 4,p<? 



a 



a 



a 



t (!) 

-?7 T (7 2 i o\ 2 

T (2) T 

V °22 °2\ 



definite body, <j^ T 



G Mat(p/2g) 



°i = °"i = °i with positive 



_ (1) (2)T 
°22 ? u 22 



-0", 



,(2) 

22 r 5 



(71 ?7 



G Mat(p/g) 

G Mat(2p/g) 
with positive definite body, a 2 = a 2 \. 



erf = <j\ with positive definite body > , 



(4.1) 
(4.2) 



o"n 0-12 V 

t T 
-if Tj 2 



a\ = (7i 



0"12 

-^12 ^11 



(4.3) 
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Also, we will use the sets 



E« g = { a e ^ pq \ 4 = a 2 ) = S« g n E Am (4.4) 
and 

^ pq ={o-e^ pq \a 2 eCXJ^(q)} (4.5) 

where CU ^' (q) is the set of the circular orthogonal (COE, (3 = 1), unitary (CUE, 
(3 = 2) or unitary-symplectic (CSE, (3 = 4) ensembles, 

AeU (2) ( 9 ) ,0 = 2 1(4.6) 

A=(r s ®l 9 )^(r/®l 9 )GU( 2 )(2g) ,0 = 4 J 

The index "f" in Eq. (14.41) refers to the self-adjointness of the supermatrices and the 
index "c" indicates the relation to the circular ensembles. We notice that the set classes 
presented above differ in the Fermion-Fermion block. In Sec. El we show that this is the 
crucial difference between both methods. Due to the nilpotence of -B's Fermion-Fermion 
block, we can change the set in this block for the Fourier-transformation. The sets of 

yO(f) 

P,pq J P,pq 



CU<« {q)= {AeG\( l2 q,C) 



matrices in the sets above with entries in A and Ai are denoted by , an d 



E £> respectively. 

The proof of the superbosonization formula [19, 20J given below is based on 
the proofs of the superbosonization formula for arbitrary superfunctions on real 
super symmetric Wishart-matrices in Ref. [TH] and for Gaussian functions on real, 
complex and quaternionic Wishart-matrices in Ref. [37] . This theorem extends the 
superbosonization formula of Ref. [20J to averages of square roots of determinants over 
unitary-symplectically invariant ensembles, i.e. /3 = 4, b = c = and d odd in Eq. (|3.7|) . 



The proof of this theorem is given in Appendix A 



Theorem 4.1 (Superbosonization formula) 

Let F be a conveniently integrable and analytic superfunction on the set of (7 2 c + 71 <i) x 
(72c + 71 d) supermatrices and 

a — c+ld — 1 /A _ S 
k = + . (4.7) 



7i 72 



With 



a > c , (4.8) 

we find 

J F(B) exp (-eStrij) d[V] = C { a ^ d J F(p)exp (-eStrp) Sdet p K d[p], (4.9) 



m. v o(c) 
where the constant is 



■ >-~\ cd VolflJ (/3) ( 
C% = (-27r Tl )- ad ( --] 2"V ac/2 * 



72/ ~ Vol (U (/3) (a-c)) 
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(7iK + 2(w-d)/, 

? 4(n-l)//3 7r 2(n-l)//3 



A 11 ,4(n-l)//3 7r 2(n-l)//3 " V*"^ 
n=l 



VFe define the measure d[V] as in Corollaru \3.£\ and the measure on the right hand side 
is d[p] = d[pi]d[p 2 ]d[ri] where 



d[pi] = ] J dp nn i X < 



n=l 



EI ^Pnml , /? = 1, 

l<n<m<c 

n dRep nm idImp nm i , (3 = 2, (4.11) 

l<n<m<c 

n ciRep nm iic/Imp nm iiciRep nm i2(iImp nm i2 , /3 = 4, 

l<n<m<c 

d[ P2 ) = FUi 4//3) |A d (e^')| 4//3 II ^f— ( 4 ' 12 ) 

n=l 

c d 

<*M = UU^mdvlJ- (4-13) 

n=l m=l 

Here, p 2 = t/diag (e^ 1 , . . . , e J</,<i ) [A, Z7 G U (d) and dfx{U) is the normalized Haar- 
measure of U (d) . We introduce the volumes of the rotation groups 

voi (u«(»)) =n^j <«■"> 

and the ratio of volumes of the group flag manifold and the permutation group 

u « -dill r(2i/,8) ' (4J5) 

T/ie absolute value of the Vandermonde determinant A^(e 1 ^) = F] (e Vn -e*" 1 ) 

l<n<m<c( 

refers to a change of sign in every single difference (e l</3n — e lipm ) with "+" if <p m < <p„ 
and wzi/i "— " otherwise. Thus, it is not an absolute value in the complex plane. 

The exponential term can also be shifted in the superfunction F. We need this 
additional term to regularize an intermediate step in the proof. 

The inequality (14. 8 p is crucial. For example, let (3 = 2 and F(p) = 1. Then, the 
left hand side of Eq. H4.9f) is not equal to zero. On the right hand side of Eq. (14.9j) . the 
dependence on the Grassmann variables only stems from the superdeterminant and we 
find 

Sdet = / det(P ' d e t 'f" T 4"l = (4.16) 

for k < d. The superdeterminant Sdet p is a polynomial of order 2c in the Grassmann 
variables {n nm ,n* m } and the integral over the remaining variables is finite for k > 0. 
Hence, it is easy to see that the right hand side of Eq. f)4.9p is zero for k < d. This 
inequality is equivalent to a < c. 

This problem was also discussed in Ref . [3T] • These authors gave a solution for the 
case that ( 14.81) is violated. This solution differs from our approach in Sec. [7J 
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5. The generalized Hubbard-Stratonovich transformation 
The following theorem is proven in a way similar to Refs. [T5l EJJ. The proof is given 



m 



Appendix B| We need the Wick-rotated set S^ 1 



,(0) 
J f3,cd 



-.0(f) 
J (3,cd' 



P,cd 



nf%%flf\ particularly 



The original extension of the Hubbard-Stratonovich transformation 



[To] |2"T| was only given for 7 2 c 
and d. 



jxd = jk. Here, we generalize it to arbitrary c 



Theorem 5.1 (Generalized Hubbard-Stratonovich transformation) 

Let F and k be the same as in Theorem \4-l\ If the inequality ft4-8\ ) holds, we have 



F(£)exp (-eStrS) d[% 



IK 



a 



acd 



F (p) exp (-eStrp) det ( e'^D^f 



S(r 2 ) 



09) 
acc( 



det p* 



r.(0) 
J /3,cd 



|A,(r 2 )|^ ^ ^ J 



|A d (e^r 2 )| 4 //3 
F(p) exp(-eStrp)| i/)=0 d[p] 



^ cd d[p] 



(5.1) 



with 



Pi 


e^/ 2 pr, 


-e#/2 p t 


^ (P2-Plpl 1 pr l ) 



(5.2) 



T/ie variables r 2 are £/ie eigenvalues of the supermatrix p 2 . T7ie measure d[p] = 
d[pi]d[p 2 ]d[r]] is defined by Eqs. (f^TTTJ) and (f-^-iffi ). For i/ie measure d[p 2 ] we take the 
definition ft4.ll ) for 4/(3. The differential operator in Eq. ft5. 1\) is an analog of the 
Sekiguchi- differential operator [38] and has the form I21\( 

1 



D 



dr2 



A d (r 2 ) 



det 



d-m 
n2 



9 / j ,2 1 
dr-. +{d - m) (3^ 



Kn,m<d 



The constant is 



a 



acd 



cd 



(2vr 7l )- ad (^ 7^° C/2 



VolfU^f 



Vol (U^(a-c))FUi 4//3) 



(5.3) 



(5.4) 



Since the diagonalization of p 2 yields an |A^(r 2 )| 4 /' 8 in the measure, the ratio of 
the Dirac-distribution with the Vandermonde-determinant is for Schwartz-functions 
on Herm (4/(3, d) well-defined. Also, the action of D^JJ 3 ^ on such a Schwartz-function 
integrated over the corresponding rotation group is finite at zero. 

The distribution in the Fermion-Fermion block in Eq. (15. lft takes for (3 G {1,2} the 
simpler form [Jj)J EJJ 

c 5(r 2 ) 



D 



dr2 



FU (4/» 

n=l 



r(a-c + l + 2(n-l)//3) 

(-7r)2(n-l)//3 r ( 7lR ) 



ll^T^T^n) • 

n=l a '«2 



(5.5) 
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This expression written as a contour integral is the superbosonization formula [39]. For 
(3 = 4, we do not find such a simplification due to the term |A(r2)| as the Jacobian in 
the eigenvalue-angle coordinates. 

6. Equivalence of and connections between the two approaches 

Above, we have argued that both expressions in Theorems 14.11 and 15.11 are equivalent 
for (3 G {1,2}. Now we address all (3 G {1,2,4}. The Theorem below is proven in 



Appendix C The proof treats all three cases in a unifying way. Properties of the 
ordinary matrix Bessel-functions are used. 

Theorem 6.1 (Equivalence of Theorems 14.11 and 15.1ft 



The superbosonization formula, \4-l\ and the generalized Hubbard-Stratonovich 



transformation, \5.1\ are equivalent for superfunctions which are Schwartz-functions and 
analytic in the fermionic eigenvalues. 

The compact integral in the Fermion-Fermion block of the superbosonization 
formula can be considered as a contour integral. In the proof of Theorem 16.11 we 
find the integral identity 

-1 p'(l-7l«)^n^ n 



J F(e^)\A d (e^)\^f[ ^ 

[0,27T] d 

tt ^-^r(i + yg) ( D (my- c F (r ) 

l = lr(2//3 + l)r( 7l K-2(n-l)//3) V dr2 J {2} 



for an analytic function F on C d with permutation invariance. Hence, we can relate 
both constants ( 14.101) and (E 



c2 L\ T (VP + i)r(7i« - 2(« - i)l P) ' ( ' ' 

The integral identity (16. ip is a reminiscent of the residue theorem. It is the analog of 
the connection between the contour integral and the differential operator in the cases 
(3 G {1,2}, see Fig. [TJ Thus, the differential Operator with the Dirac-distribution in 
the generalized Hubbard-Stratonovich transformation restricts the non-compact integral 
in the Fermion-Fermion block to the point zero and its neighborhood. Therefore 
it is equivalent to a compact Fermion-Fermion block integral as appearing in the 
superbosonization formula. 

7. The general case for arbitrary positive integers a, b, c, d and arbitrary 
Dyson-index (3 G {1,2,4} 

We consider an application of our results. The inequality (14 ,8p reads 

N > 7l A; (7.1) 
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lm 




Figure 1. In the superbosonization formula, the integration of the fermionic 
eigenvalues is along the unit circle in the complex plane (dotted circle). The eigenvalue 
integrals in the generalized Hubbard-Stratonovich transformation are integrations over 
the real axis (bold line) or on the Wick-rotated real axis (thin line at angle t/j) if the 
differential operator acts on the supcrfunction or on the Dirac-distribution at zero 
(bold dot, 0), respectively. 



for the calculation of the /c-point correlation function (12.41) with help of the matrix Green 
function. For (3 = 1 , a N x iV real symmetric matrix has in the absence of degeneracies 
N different eigenvalues. However, we can only calculate fc-point correlation functions 
with k < N/2. For N — ► oo, this restriction does not matter. But for exact finite N 
calculations, we have to modify the line of reasoning. 
We construct the symmetry operator 
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(7.2) 



from (mi + m-i) x {n\ + ri2) supermatrix to (rri2 + mi) x (ri2 + ni) supermatrix. This 
operator has the properties 

(7.3) 
(7.4) 
(7.5) 



(7.6) 



Let a, b, c, d be arbitrary positive integers and (3 G {1,2,4}. Then, the equation 
(17.6p reads for a matrix product of a (72c + 71 d) x (0 + 716) supermatrix with a 
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(0 + jib) x (72c + jid) supermatrix 



St 



& 



St 



6 



6 



St 



(7.7) 



With help of the operator (5, we split the supersymmetric Wishart-matrix B into two 
parts, 

B = Bi + S(B 2 ) (7.8) 

such that 



t 



(7.9) 



The supervectors 6 f^j?) are °f the same form as ■ Let a be a quadratic 
supermatrix, i.e. mi = m and m 2 = n 2 . Then, we find the additional property 

(-l) m2 Sdet~V. (7.10) 



Let S 



(0) 
I3,pq 



Sdet6(a) 



-0(c) 



6 £ 



-0(c) 



and the Wick-rotated set 



n^E^°p (? n^ C ' ) . Then, we construct the analog of the superbosonization formula and 
the generalized Hubbard-Stratonovich transformation. 

Theorem 7.1 



Let F be the superfunction as in Theorem \4-l\ and 
a — c + 1 b — d + 1 



K 



71 



72 



and 



b = b + j 2 e 



(7.11) 
(7.12) 



y4/so ; Zei e G N and 

a = a + 7ie 

a > c 6 > d. (7.13) 

We choose the Wick-rotation e 1 ^ such that all integrals are convergent. Then, we have 

F(^)exp (-eStrB^) d[V) = 



7i 



72 ec 



72 



71 eci 



F(5 v ,)exp (-eStrBj d[V] 



a 



SF 



d[p (2) ]d[p (1) ]F(p( 1 ) + e"V 2) ) exp [-eStr (p (1) + e^p (2) )] x 



„0(c) SO(c) 
^ a — 1 a 



^0,cd ^4/f3,dc 

xSdet ^W^Sdet^W 2 ) 



C 



11s 



d[p^]d[p^ 



5 ( r« 



(7.14) 



E (0) g(0) 



/3,cd ^4/f3,cd 



4//3 



det" +6/72 pf ) det o/7l - K ^ 2) x 



P2 
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x D 



)(4//3) \ 



D 



(2) 



b-d 



F(pW + e^'p {2) ) exp [-eStr (p (1) + e^p 



(2) 



where the constants are 

C Hg = r_l\d(a-c) e itp(ad-bc) 

Here, we define the supermatrix 



2 \ 72ec 

7i 

2 \ 72ec 

71/ 



71 ed 



r (J3) r W0) 



72/ 



71 ed 



acd Icic ■ 









P 2 1} - rf'VT W + ^pf . 



(7.15) 

(7.16) 
(7.17) 

(7.18) 



The set 91 is given as in corollary \3.2l The measures d[p^>\ = d[p^]d[p^]d[rj\ and 



d[p^} = d[pi]d\fi2 ]d[fj] are given by Theorem \4-l\ The measures ft4.ll]) for (3 and 4/ '(3 
assign d[p± ] and d[p^] in Eqs. ( |7. O ) and ( |7. 15]) , respectively. In Eq. ft7.14\ ), d[p^] 
and d[p^\ are defined by the measure U-IB) for the cases (3 and 4//3, respectively, 
and, in Eq. ft7.15]) , they are defined by the measure ft4-H\ ) for the cases 4/ '(3 and (3, 
respectively. The measures d[rf\ and d[fj] are the product of all complex Grassmann pairs 
as in Eq. ( gTTgp . 

Since this Theorem is a consequence of corollary 13.21 and Theorems 14.11 and 15. 1[ the 

proof is quite simple. 

Proof: 

Let e G No as in Eq. (17.121) . Then, we use corollary 13.21 to extend the integral over V to 
an integral over V. We split the supersymmetric Wishart-matrix B as in Eq. ( 17. 81) . Both 
Wishart-matrices Bi and B 2 fulfill the requirement H4.8|) according to their dimension. 
Thus, we singly apply both Theorems 14. 1 1 and 15 . 1 1 to B\ and B 2 . □ 
Our approach of a violation of inequality (14.81) is quite different from the solution 
given in Ref. [5T]. These authors introduce a matrix which projects the Boson-Boson 
block and the bosonic side of the off-diagonal blocks onto a space of the smaller 
dimension a. Then, they integrate over all of such orthogonal projectors. This integral 
becomes more difficult due to an additional measure on a curved, compact space. We 
use a second symmetric supermatrix. Hence, we have up to the dimensions of the 
supermatrices a symmetry between both supermatrices produced by &. There is no 
additional complication for the integration, since the measures of both supermatrices 
are of the same kind. Moreover, our approach extends the results to the case of (3 = 4 
and odd b which is not considered in Ref. [51] . 



8. Remarks and conclusions 



We proved the equivalence of the generalized Hubbard-Stratonovich transformation 
[15. 21J and the superbosonization formula [19. 20j. Thereby, we generalized both 
approaches. The superbosonization formula was proven in a new way and is now 
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extended to odd dimensional supersymmetric Wishart-matrices in the Fermion-Fermion 
block for the quaternionic case. The generalized Hubbard-Stratonovich transformation 
was here extended to arbitrary dimensional supersymmetric Wishart-matrices which 
not only stem of averages over the matrix Green functions. [HI [2H EH 122] Furthermore, 
we got an integral identity beyond the restriction of the matrix dimension, see Eq. ( 14. 8ft . 
This approach distinguishes from the method presented in Ref. [31 J by the integration 
of an additional matrix. It is, also, applicable on the artificial example (3 = 4 and odd 
b which has not been considered in Ref. [3l] . 

The generalized Hubbard-Stratonovich transformation and the superbosonization 
formula reduce in the absence of Grassmann variables to the ordinary integral identity 
for ordinary Wishart-matrices. [2H1 [20] In the general case with the restriction (14. 8p . 
both approaches differ in the Fermion-Fermion block integration. Due to the Dirac- 
distribution and the differential operator, the integration over the non-compact domain 
in the generalized Hubbard-Stratonovich transformation is equal with help of the residue 
theorem to a contour integral. This contour integral is equivalent to the integration over 
the compact domain in the superbosonization formula. Hence, we found an integral 
identity between a compact integral and a differentiated Dirac-distribution. 
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Appendix A. Proof of Theorem 14.11 (Superbosonization formula) 

First, we consider two particular cases. Let d = and a > c be an arbitrary positive 
integer. Then, we find 

B e S° jC0 = = Eg C Herm c). (A.l) 

We introduce a Fourier-transformation 

J F(B) exp (-etrlj) d[V] = 

= (g)' (2^ C(C " 1)/2 J J jr F{ai) exp (i %T Bot) diVWr] (A.2) 

Herm (/3,c) U0 ac 

where the measure d[o~\\ is defined as in Eq. (14. lOf) and a± = o\ + zel 72C . The Fourier- 
transform is 

TF{a x )= J F(pi)exp(-ztrp 1 (7 1 )d[ / 9 1 ]. (A.3) 

Herm (/3,c) 
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The integration over the supervectors, which are in this particular case ordinary vectors, 
yields 

/ + \ — a /7i 

exp (i tr Bat) = det f — J • ( A - 4 ) 



R /3ac 

The Fourier-transform of this determinant is an Ingham-Siegel integral [lOl W\\ 

J exp(-ztrpio-i)det (-2o- 1 + ) _a/7l rf[o-i] = Gf2 cc det p* exp (-etrpx) 0(pi), (A.5) 

Herm (/3,c) 

where the constant is 

<*--(?r n S (A - 6) 

j=a-c+l y 



and the exponent is 



k = . (A.7) 

7i 72 

T(.) is Euler's gamma-function. This integral was recently used in random matrix theory 
[29] and is normalized in our notation as in Ref. [21]. Thus, we find for Eq. (1A.2I) 

J F(S)exp(-etrs)d[^ = c2 J F{p) exp (-etrpi) detp?d[pi], (A.8) 



r 0(c) 
/3,c0 



which verifies this theorem. The product in the constant 

Vol fu W (a)) 
Vol fU (/3) (a-c)J 

is a ratio of group volumes. 

In the next case, we consider c = and arbitrary d. We see that 

SeSjL (A.10) 
is true. We integrate over 



F(B)exp(etTB)d[V], (A.ll) 



A 



2nd 



where the function F is analytic. As in Ref. [19], we expand F(B) exp (etx B^j in the 

entries of B and, then, integrate over every single term of this expansion. Every term 
is a product of 5's entries and can be generated by differentiation of (tr ABj with 

respect to A 6 g d for certain nGN. Thus, it is sufficient to proof the integral theorem 
J (tr ABj d[V] = C { 2 J (tv Ap 2 ) n det p^d[p 2 }. (A.12) 
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Since Sjgod is generated of by analytic continuation in the eigenvalues, it is 

convenient that A G E°$ d . Then, A~ 1/2 is well-defined and A~ l / 2 p 2 A- 1 / 2 G S^- 
We transform in Eq. ( 1A.13I) 

v _> a- 1 ^!? ; f>t _ and P2 ^ A- 1 / 2 p 2 ^" 1/2 - 

The measures turns under this change into 

d[V] -> detA a/71 d[V] and 
d[p 2 ] -> detA- K+a ^M[p 2 ], 

where the exponent is 

a + 1 d- 1 

7i 72 

Hence, we have to calculate the remaining constant defined by 

,03) 



trfl d[V] = CZ / (trp 2 )"detp 2 -^[P: 



.-0(c) 



(A.13) 

(A.14) 
(A.15) 

(A.16) 
(A.17) 



This equation holds for arbitrary n. Then, this must also be valid for F(B) = e = 1 in 
Eq. flA~TT]) . The right hand side of Eq. flATH is 



exp (trS) d[V] = (-2tt) 



(A.18) 



On the left hand side, we first integrate over the group U ( - 4///3 - ) (d) and get 
exp (trp 2 ) detp2 K d[p 2 ] = 



-0(c) 



FUj //3) / |A d (e^)| 4//3 n ex P^ e ^") e " 

[0,27T] d 



2tt ' 



(A.19) 



(A.20) 



We derive this integral with help of Selberg's integral formula [14] . We assume that 
— A/f3 and 71 k are arbitrary positive integers and (3 is even. Then, we omit the 
absolute value and Eq. (IA.20I) becomes 

exp(trp 2 )detp 2 - K rf[p 2 ] = FUf Aj ^ ' 



-0(c) 
J /3,(M 



tiSaJAI r (71 /«) 



(A.212 



A=l 



We consider another integral which is the Laguerre version of Selberg's integral [H] 



^d( X ) II eX P (-7l^n) X^dXn = A 



n=l 



1 9 



7i 



n r (^+i) (71 a, 



-€-1 



'3 / n=l 



A=l 



J_ r(l + n/3/2)r(e + l + (n-l)/3/2 



n=l 



7i 



+l+/3(d-l)/2 1 



r ( 1 + 0/2 



(A.22) 
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where £ is an arbitrary positive integer. Since (3 is even the minus sign in the 
Vandermonde determinant vanishes. The equations (IA.21I) and ()A.22j) are up to the 
Gamma-functions polynomials in k and £. We remind that ( 1A.22I) is true for every 
complex £. Let Re£ > 0, we have 



]__d_ 

7i d\ 3 



n 



-t-i 



7i 



•d(e+l+/3(d-l)/2) 



ir^+l + (n-l)/3/2 
r f 1 + n/3/2 



A=l 



n 



< const. < oo and 



< const. < oo. 



(A.23) 



(A.24) 

"ir(£ + i)r + ' 

The functions are bounded and regular for Re£ > and we can apply Carlson's theorem 
[H]. We identify £ = — j%k and find 

exp (trp 2 ) detp2 K d[p 2 ] = 



.-0(c) 



rf tt/m A r 1 + n/? / 2 r 1 - 7i k + (n - 10/2 
= 7i ad FU? //3) J] ^ ^. (A.25) 

n=i r(i + /3/2jr( 7 iK)r(i-7iK) 

Due to Euler's reflection formula T(z)T(l — z) = it/ sin(irz), this equation simplifies to 



exp (trp 2 ) det p 2 K d[p 2 ] = ^FV^^ J] — - 



aWVfl A ^ 4(n - 1)//3 r (1 + 2n/(3) 



v 0(c) 
i8,0d 

or equivalent 



2/9<i(d-l)/2 



[0,2tt 
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d l<n<m<(i 



sin 
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P d 



2//3)r( 7l «-2(n-l)//3) 



(A.26) 



JJexp ( lie tlpn ) e~ llpna 



n=l 



difn 

2vr 



r 1 + n/3/2 



n=i T 1 + /3/2 r o + 1 + (n - l)/3/2 



(A.27) 



Since a is a positive integer for all positive and even /3, the equations above are true for 
all such /3. For constant natural numbers a, d, 71 and complex (3 with Re/3 > 0, the 
inequalities 



0,2tt] 



< 



[0,27r] d 
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d l<n<m<e£ 
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Kn<m<d 



sin 



sin 



P d 

J^exp ( 7l e^)e-^ 

n=l 
Re/3 rf 

exp (71 cos w 

n=l 



r 1 + n/3/2 



1 T 1 + p/2 r a + 1 + (n - 1)0/2 



2tt 

2tt 

< 



< 



< 00 and 



(A.28) 
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< const. 2 - Rc ^ d - 1)/2 < oo (A.29) 
are valid and allow us with Carlson's theorem to extend Eq. ()A.27j) to every complex (5, 



in particular to /3 — 1. Thus, we find for the constant in Eq. (1A.17I) 

. II;/- I ).. )' _2i // — I )/ >' n ' 

, —ad 



a 



aOd 



-27T7 



n 



t 4(n-l)/f3 n 2(n-l)/f3 



llr(a + l + 2(n-l)//3) 
Now, we consider arbitrary d and a > c and split 



~3 

Xj 



(A.30) 



(A.31) 



and 



1 3=1 



a f 



E 

i=i 
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a f 1 
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3=1 
(C) 
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7 



B 2 \ B 22 



(A.32) 



such that Xj contains all commuting variables of and Xj depends on all Grassmann 
variables. Then, we replace the sub-matrices .612,-621 and -B22 by Dirac-distributions 

F(B)exp f-eStrB) d[V] = 



Ci 



d[r))d[fj]d[V)d[p 2 ]d[a 2 }F 



Herm (4//3,d)2 £R (A 2cd )2 



X 



xexp 
where 



eStr B - % (tr(pj + S 21 )^ + tr <tJ(p„ - S 12 ) - tr(p~ 2 - #22)^2 



2vr\ 
7 / 



2cd 



7l s 2d(d-l)//3 / 7l 
W V27T 



(A.33) 



(A.34) 



The matrices p,, and 0^ are rectangular matrices depending on Grassmann variables 
as in the Boson-Fermion and Fermion-Boson block in the sets (14. 1114.31) . Shifting 

xj Xj + {°zy 4 x i and x) -> x] - (<4)~ > we § et 



F(5)exp (-eStrS) 



dMrf[r/]d[y]d[p2]d[a 2 ]F 



Herm(4//3,d)2 <K (A 2cd )2 



x 



xexp 



eStr B — 1 ( tr atBna^ (a 2 ) 1 + tr p]^ + tr crip,, - tr(p 2 - B 2 2)o- 2 



(A.35) 
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This integral only depends on Bu and Bw Thus, we apply the first case of this proof 
and replace Bu. We find 



F(B)exp (-eStvB) d[V) = 

d2>C x f f f dWWW[ P i]d\p 2 ]d[a 2 ]F 

Herm(4//3,d) 2 9t (A 2cd ) 2 



B\\ p v 

~P\ P2 



det pi x 



xexp 



e{tiB 22 - trpx) + » (tro-tpicr^ (ovf) 1 - trpjo^ - tr o\p n + tr(p 2 - #22)^2 J A.36) 



with the exponent 

a — c + 1 1 . t , 

k = . (A.37 

7i 72 

After another shifting — ► cr^ — p\ l p^t and cr| — » cr| — a 2 f p\Pi 1 '■> we integrate over 
d[f/] and I? 22 and have 

exp f-eStr B] d[V] = 



d[rj\d[pi}d[p 2 }d[(T2}F 



yO(c) Herm (4//3,d) 2 A 2cd 
/3,c0 

+ n (a-c)/7i 



Pi Pr? 
-Pj P2 



x det p\ det (a^ ) 71 exp [— e tr p x + 1 (tr pj ? p 1 1 p r? cr^ + tr p 2 a 2 )] 

where the exponent is 

a — c+1 d — 1 



K 



7i 

and the new constant is 
C 2 ~- 



+ 



72 



(A.38) 
(A.39) 



l \ ad /27T\ Cd / 7l \ 2d(d-l)//3 / 7l \d 

2tt/ VT^y W/ V2tt/ 

We express the determinant in as in Sec. [2] as Gaussian integrals and define a 
new (72(0 — c) + 0) x (0 + 71 d) rectangular supermatrix V„ ew and its corresponding 
(0+jid) x (O + 71G?) supermatrix B new = 7 _1 KewK]ew Integrating <r 2 and p 2 , Eq. (1A.38I) 
becomes 



(A.40) 



exp -eStr B = 7 " cd C ( 



~-cdri(P) 

J acQ 



S 0(c) A 2(a-c)d 



Pi 


Pv 


-4 


-B new PrjPl Prj 



exp ( -etrpi +etr(B ncw - ?7 f pi S) ) det Pid[V new ]d[rj\d\pi]. 



(A.41) 



Now, we apply the second case in this proof and shift p 2 — > P2 + plpi Pr; by analytic 
continuation. We get the final result 



-F(-B) exp ( -eStr B ) d[V] = Cf-j / F (p) exp (-eStr p) Sdet p K d[p] 



(A.42) 



.-.0(c) 
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with 

^acd — I °acO°a-c,Od 



-27T 7 l) 



cd 



2tt 

72 



Vol (U W(a)) - r( 7l « + 2(n-d)/)9) 



-2d(d-l)//3 \ Z7T ) I 

(_2)(c-a)rf2 c 



Vol (U ^(a-c 
2 d2 Vol (lJ (1) (a 



n=l 



Vol (lJ {1) (a-c + 2d) 
Vol(U (2) (a) 



Vol (U {2) (a-c + rf) 
2 -(2a+i-c) C y ol (\j (1) (2a + 1) 

Vol (U (1) (2(a-c) 



^4(n-l)//3 7r 2(n-l)//3 
, /3 = 1 

, /9 = 2 



= 4 



(A.43) 



Appendix B. Proof of Theorem 15.11 (Generalized Hubbard-Stratonovich 
transformation) 

We choose a Wick-rotation e*^ that all calculations below are well defined. Then, we 
perform a Fourier transformation 



F(B) exp ( -eStr B I d[V] = d / / .FF(cr) exp ( zStr 5(7+ ) d[V]d[<r 



where cr + = a + ^£ll 



72C+7ia!) 



FF(a) = / exp (-zStr pa) d[p] , 



and the constant is 



Ci 



2vr \ 2cd / 72 \ c / 72 \/ 3c ( c - 1 )/ 2 / 7 i\ d / 7l \2d(rf-i)//3 



7 / \27T/ V 7T 

The integration over yields 



27T/ V 7T 



exp ( -eStr B ) d[V] = C 2 I J r F( ( x)Sdet ~ ahl a + d[a] (B.4) 



with 



7 J \2tT/ V 7T / V27T/ V 7T / V27T 

We transform this result back by a Fourier-transformation 

exp (-eStr b) = C7 2 / F(p)/^' a) (p) exp (-eStrp) 

^f3,cd 



(B.l) 



(B.2) 



(B-3) 



(B.5) 



(B.6) 
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where we have to calculate the supersymmetric Ingham-Siegel integral 

I^ a \p)= J exp (-*Strpa+) Sdet ~ ahl a + d[a}. (B.7) 



v(-vo 



This distribution is rotation invariant under U^\c/d). The ordinary version, d = 0, of 
Eq. dHSD is Eq. flA3|) . 

After performing four shifts 



eri -»• eri - 0jj (02 + «e # £l 7l( i) 1 crt, 
o-fj ^ o-fj- p{ x Pn (°"2 + ze # el 7ld ) , 
4 - 4 - ((7-2 + ie^et ltd ) p\p{\ 



and defining 



we find 



P2 -> P2 - pJPi Vr,, 



Pi 



e^Pr? 



* (pa - pJprVJ 



(B.8) 
(B.9) 
(B.10) 
(B.H) 



(B.12) 



F(B)exp (-eStrfl) d[V] = C 2 / F (p) I{p) exp (— eStr p) d[p], (B.13) 



where 
?(p) = 



exp 



eStrp - 2 (trpi0! - trp 2 2 + tr 0~pi 0^(02 + «e^el 7ld ' 



-1 



x 



'.(->/>) 



X 



det(e-^0 2 + ?£l 7ld ) 



a/71 



<f[4 



(B.14) 



det(ai + 2£l 72 c) 

We integrate over d[fj] and apply Eq. (IA.5I) for the G?[0i]-integration. Then, Eq. flB.141) 

reads 

7(p) = C i 3e~^detp^e(p 1 )x 

x J exp(-«trp 2 (02 + ^U 7ld ))det(e-^0 2 + ^) (a " c)/7l rfN (B.15) 



Herm (4//3,<f) 

with the constant 



3 = * 



-/3ac/2 / 7 

2™ 



cd 



G, 



(B.16) 



see Eq. ( 1A.6I) . The exponent k is the same as in Eq. (14.71) . As in Ref. [21], we decompose 
02 in angles and eigenvalues and integrate over the angles. Thus, we get the ordinary 
matrix Bessel-function 



exp (itr r 2 Us 2 U^) dp(U) 



(B.17) 



U< 4 /«(d) 
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in Eq. (IB. 151) which are only for certain (3 and d explicitly known. However, the analog 
of the Sekiguchi differential operator for the ordinary matrix Bessel-functions DfJ^\ 
see Eq. (15.31) . fulfills the eigenvalue equation 

(B.18) 



D 



dr2 



^ /P, (r 2 ,s 2 ) = ^ i rdets 1 2 ^^\r 2 ,s 2 ). 



Since the determinant of a 2 stands in the numerator, we shift a 2 — > o 2 — ze l %l 71( 2 and 

* A ov on inf orrrnf mn fxt rov iT„ 

dr-2 



replace the determinants in Eq. ( IB. 151) by D^J^ . After an integration over a 2 , we have 



I( P ) = C,e-^det P «e( Pl ) (e-^D^j 



The constant is 



2d(d-l)/f3 / 27r \ °' 



\7i/ FU 



(4//3) ' 



Summarizing the constants (1B.5[) and ( IB. 201) . we get 

«i Vol(U (/3) (a 



^2 = W4 = 2- c (2vr 7l )- ad (^) 7^ 



Vol(U^(a-c))FUi 4//3) 



(B.19) 



(B.20) 



(B.21) 



Due to the Dirac-distribution, we shift D^JJ 3 ^ from the Dirac-distribution to the 
superfunction and remove the Wick-rotation. Hence, we find the result of the Theorem. 



Appendix C. Proof of Theorem 16.11 (Equivalence of both approaches) 

We define the function 



Pi 



-P\ 



Ur 2 W - plpV 1 



PvPi Pv 



X 



FW = J J f 

U 4 //3(rf) Herm (/3,c) A 2cd 

x exp [— e(trpi - tr(r 2 - P^prV??)] detVi^MPi^M^)- (C.l) 
Then, we have to prove 



<* e l(l—K,)tp n 



[0,2tt] 



n=l 



2tt 



• (C2) 

r 2 =0 



Since F is permutation invariant and a Schwartz-function, we express F as an integral 
over ordinary matrix Bessel-functions, 



F(r 2 )= / g( q )^ /P \r 2 ,q)\A d (q)\^dq, 



(C.3) 



where g is a Schwartz-function. The integral and the differential operator in Eq. (10.21) 
commute with the integral in Eq. (1C.3I) . Thus, we only need to prove 

v r)(e-,,)iA J( e-)i^ri- 1 " 7 '*"^ 

n=l 



a 



09) 



[0,27r] d 

'■0 (f_i\*r>WflY 



2tt 



(C.4) 



r 2 =0 
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for all q G Sf where S\ is the unit -circle in the complex plane. The right hand side of 
this equation is with help of Eq. (IB. 181) 

D% m )*~ C <P? fn {r 2 ,q) = (-z 7 i) d(a - c) detg(-)^. (C.5) 

7 r 2 =0 

The components of q are complex phase factors. The integral representation of the 
ordinary matrix Bessel-functions (IB. 171) and the d[(p] -integral in Eq. (10.41) form the 
integral over the circular ensembles CU (d) . Thus, q can be absorbed by e' l ' Pj and 
we find 



J ^(e^,q)\A d (e^H ^ 

[0,27r] d n=1 

/ j „l(l— 'YlK.)ip n J , „ 
^/P) (eWf 1} | Ad (eW) |V/» JJ £ (c . 6) 



[0,27r] d n_1 



The ordinary matrix Bessel-function is at q = 1 the exponential function 



(ew.^ap k^ . (C.7) 



71=1 



With Eq. ( 1 A. 271) . the integral on the left hand side in Eq. ( 1C.6I) yields with this 
exponential function 

d ^(l-Ti^iPnexp ( !Tl e^) rf^ n 



/ |A d (e^)i 4//3 Il 



[0,2tt]' 



2tt 

n=l 



= ^ TT ^-^r(i + 2n//?) 

1 n; llr(l + 2//3)r(a-c+l + 2(n-l)//5) 

= pu (4//J) 11 r(a-c+l + 2(n-l)//3)- (Ci 
Hence, the normalization on both sides in Eq. (1C.2I) is equal. 
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